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Abstract 

Based on the low energy effective Hamiltonian with naive factorization, we calculate the 
branching ratios(BRs) and CP asymmetries (CPAs) for the twenty three double charm decays 
B/B s — )■ Z?|*jl)|*j in both the standard model (SM) and the minimal supergravity (mSUGRA) 
model. Within the considered parameter space, we find that (a) the theoretical predictions 
for the BRs, CPAs and the polarization fractions in the SM and the mSUGRA model are 
all consistent with the currently available data within ±2cr errors; (b) For all the considered 
decays, the supersymmetric contributions in the mSUGRA model are very small, less than 7% 
numerically. It may be difficult to observe so small SUSY contributions even at LHC. 

PACS numbers: 13.25.Hw, 14.40.Lb, 12.60.Jv, 12.15.Ji 
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I. INTRODUCTION 



Within the standard model (SM), the double charm decays of B u d and B s Mesons 
considered here are dominated by the color-favored "Tree" transition h — > ccd(s), while the 
color-suppressed "Penguin" transition is generally small. If the penguin contribution was 
absent, the mixing induced CP asymmetry (CPA), denoted as 5*/, would be proportional 
to sin(2/3), while the direct CPA, denoted as C/, would be zero. In some new physics 
models beyond the SM, the penguin contributions can be large and may change the SM 
predictions for the branching ratios and the CP asymmetries (CPA) significantly. The 
study of these double charm B/B s meson decays therefore plays an important role in 
testing the SM as well as searching for the signals of the new physics (NP). 

Experimentally, the BaBar and Belle Collaboration have reported the measurement of 
the direct CPA in B° — > D + D~ decay 

C(f}° — ^ n + D~} f -0-91 ±0.23 ±0.06 (Belief), 

~^ u u ) — | -0.07 ±0.23 ±0.03 (BaBar [2]). [L) 

It is easy to see that Belle found an evidence of CP violation in B° — > D + D~ at the 
4.1a level [l[, but BaBar did not [2]. On the other hand, such a large direct CPA in 
B° — > D + D~ decay has not been observed in the measurements for other similar decay 
modes: such as B° -> D M+ D w ~ , B~ D Wa D M ~ and B° s -> D ( * )+ D M ~ although 
they have the same flavor structures as B° — > D + D~ at the quark level. In the SM, the 
direct CPA's should be naturally very small in size because the penguin contributions are 
small. If the large CP violation in B° — > D + D~ from Belle is true, it would establish the 
presence of new physics. 

Up to now, by using the low-energy effective hamiltonian and various factorization 
hypothesis, many investigations on the decays of B to double-charm states have been 
carried out in the framework of the SM [Io|, ll|| or some popular new physics models 

In this paper, we will present our systematic calculation of the branching ratios and 
CP violations for double charm decays B/B s — > D^)D^} in the minimal supergravity 
(mSUGRA) model [16]. In the framework of the mSUGRA model, the new physics 
contributions to the semileptonic, leptonic and radiative rare B decays and the charmless 
two-body B- meson decays have been investigated in previous works [l7-21|. For the 



two-body B — > M\M.2 decays, the new physics part of the Wilson coefficients Ck(k 
3, ■ • • , 6),C*7 7 and C$ g in the mSUGRA model can be found in Ref. 21 



The usual route to calculate the decay amplitude for non-leptonic two-body B decays 
is to start from the low energy effective Hamiltonian for AB = 1 decays. With the 
operator product expansion method, the relevant AB = 1 effective Hamiltonian can be 
factorized into the Wilson coefficients Ci(fi) times the four-quark operators Qi(f^). As to 
Cj(/x), they have been evaluated to next-to-leading order with the perturbation theory 
and renormalization group method. The remanent and also intractable problem is to 
calculate the hadronic matrix elements of these four-quark operators. Up to now, many 
methods have been put forward to settle this problem, such as the naive or generalized 



factorization approach [22|, |23j, QCD factorization approach (QCDF) [24j . |25| | and the 
perturbative QCD (PQCD) approach (26[. For the strong phase, which is important for 
the CP violation prediction, is quite sensitive to these various approaches, and different 
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approaches may lead to quite different results. In this paper, we will use the naive 
factorization method, which is expected to be reliable for the color-allowed amplitudes, 
which are dominant contributions in these double charm decays. 

This paper is organized as follows. In the next section we will give a brief review 
for the mSUGRA model. In Sec. 1111} we introduce the basic formulas for calculating 
the branching ratios, the polarization fractions and the CP violation in the considered 
B/B s — > £)/*jZ)/*j decays. In Sec. HV] we present the numerical results for the double 
charm decays of B-meson in both the SM and the mSUGRA model. The conclusions are 
included in the final section. 



II. OUTLINE OF THE MSUGRA MODEL 



In the minimal supersymmetry model (MSSM), the most general superpotential takes 
the form 0, Ht) 



fuijQiHl^Uj + foijH^Q^Dj + fEijH^L^Ej - fiH^H^ 



(2) 



a set of terms which explicitly but softly break SUSY should be added to the supersym- 
metric Lagrangian. A general form of the soft SUSY-breaking terms is given as 

- C soft = (rn 2 Q ).. q^qLj + (ml).. u* m u Rj + (m 2 D ).. d* Ri d Rj + (m 2 L ) .. t^ Lj 

+ (™!)y e* m e Rj + A^h+h + A\h+h 2 

Amj&^u* + A Dij h?f Li d* Rj + A Eij h<?P Li e* + BuhZhg 



+^m B BB + ^m w WW + ^m d GG + H.C. 



(3) 



where qu, u* Ri , d* Ri , In, e* Ri , hi and h 2 are scalar components of chiral superfields Qi, Ui, 
Di, Li, Ei, Hi and H 2 respectively, and B, W and G are U{1)y, SU{2)l, and SU(3)c 
gauge fermions. 

In order to avoid severe phenomenological problems, such as large flavor changing neu- 
tral currents (FCNC), unacceptable amount of additional CP violation and so on, a set of 
assumptions are added to the unconstrained MSSM in the mSUGRA model. One under- 
lying assumption is that SUSY-breaking occurs in a hidden sector which communicates 
with the visible sector only through gravitational interactions. The free parameters in the 
MSSM are assumed to obey a set of boundary conditions at the Grand Unification scale 
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(4) 



where = g 2 /(4ir), while g^ (i=l,2,3) denotes the coupling constant of the U(l)y, 
SU{2) L , SU(3)c gauge group, respectively. Besides the three parameters mi, m and 
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TABLE I: Two typical sets of SUSY parameters to be used in the numerical calculation. 



CASE 


mo mi 

2 


A 


tan j3 Sign[fi] 


R 7 


A 


300 300 





2 


1.10 


B 


369 150 


-400 


40 + 


-0.93 



Aq, the bilinear coupling B and the supersymmetric Higgs(ino) mass parameter fi in the 
supersymmetric sector should also be determined. By requiring the radiative electroweak 
symmetry-breaking (EWSB) takes place at the low energy scale, both of them are ob- 
tained except for the sign of /i. At this stage, only four continuous free parameters and 
an unknown sign are left in the mSUGRA model 

tan/3, mi , m , A , sign(n). (5) 
According to the previous studies about the constraints on the parameter space of the 



mSUGRA mo del |2ll |2 814321]. we choose two sets of typical mSUGRA points as listed in 
Table E 



III. EFFECTIVE HAMILTON! AN AND OBSERVABLES 

In this section, we will give a brief review of the theoretical framework of the low energy 
effective Hamiltonian and the factorized matrix elements as well as the decay amplitudes 
for AB = 1 decays. 



A. Effective Hamiltonian in the SM and mSUGRA model 

In the SM, the low energy effective Hamiltonian for AB = 1 transition at a scale fi is 
given by [33| 
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^eff 1 — ~tk ^iQi + C2Q2 + CiQi + C 7l Q 7l + CggQ8g />+/i. c. 

p=u,c i=3 ' 



(6) 



here X p = V p bV* q for b — > q transition (p e {u,c},q G {d,s}). The detailed definition of 
the operators can be found in Ref. [33]. Within the SM and at the scale Mw, the Wilson 
coefficients Ci(Mw), • ■ ■, Ciq(Mw), C 7l {Mw) and Cs g (Mw) have been given, for example, 



in Ref. [33J. By using QCD renormalization group equations, it is straightforward to 
run Wilson coefficients Ci{M\y) from the scale /1 ~ 0(My/) down to the lower scale 
fi ~ 0{m b ). 

In the mSUGRA model, there are four kinds of SUSY contributions to the b — > d(s) 
transition at the one-loop level, depending on the virtual particles running in the penguin 
diagrams: 

(i) the charged Higgs boson H ± and up-type quarks u, c, t; 



4 



(ii) the charginos Xi 2 an d the up-type squarks u, c, t; 

(iii) the neutralinos Xi 2 3 4 an d the down-type quarks d, s, 6; 

(iv) the gluinos g and the down-type quarks d, s, b. 

In general, the Wilson coefficients after the inclusion of various contributions can be 
expressed as 

Cifrw) = C? M + Gf~ + Cf + cf + Cf, (7) 

where Cf , Cf , and Cf denote the Wilson coefficients induced by the penguin di- 
agrams with the exchanges of the charged Higgs H^, the chargino Xi2> the rieutralino 
X° 2,3,4 an d the gluino q, respectively. The detailed expressions of these Wilson coefficients 
can be found in Ref. 12 ll. 



B. Decay amplitudes in naive factorization 



The decay amplitudes of B — > D ( * ] ' in the SM within the naive factorization can 
be written as [22I 



(b -+ d w d?) = ^ ( \ c ai k + < + a< + <)}) ^ 

\ p=u,c / 



where the coefficients af = ^Cj + ^7 J + Pf with the upper (lower) sign applied when i 

is odd (even), and Pf account for penguin contributions. The factorization parameter £ 
in Eq. OH]) arises from the transformation of (V — A) (V + A) currents into (V — A) (V — A) 
ones for the penguin operators. It depends on properties of the final-state mesons involved 
and is defined as 



+ 



2ra- 



(rn c +in q )(rn b -m c ) 





2m- 



(m c +m q )(m b +m c ) 





(DD q ), 
(DD* q ), 

(D*D*). 



(9) 



BD (*) D (*)j in Eq. ([8]) is the factorized matrix element. For B — >■ D D q 



The term A 
decay mode, it can be written as 



A [bd^\d^] = ( D l W(l-7 5 )c|0) {D""\c^(l- l5 )b\B 



The decay constants and form factors [22], |34| are usually defined as 

(D q (p Dg )\qYl 5 c\0) = ~if DqP ^ 
(D* q (p Di )\q^c\0) = f j? 



(D(p D )\c^b\B(p B )) 



m B — m. 



m B — m 



(10) 

(11) 
(12) 

Fi(q 2 ), (13) 
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(D*(p D *,e*)\c^b\B(p B )) = J nq 2 we^, 



+ m 



(14) 



(L>*( Pl5 ,, £ *)|c7 M 75fc|S(p fl )) 



e* (m B + m D . )A X (q 2 ) - (p B + p D , ) fi {e* ■ p B ) 



m B + m n „ 



2m, 



~[A 3 (q 2 ) - A (q 2 )}, 



(15) 



where q = Pb — P D (*)- in terms of decay constants and form factors, the matrix element 
A IDn (,) „(*), can be written as follows 

[BU ,U q \ 



' ifD q (ml-m 2 D )F Q (m 2 Dq 
2fD*mB\Pc\Fi(m 2 

-2fD q m B \p c \A(m 2 D ), 



A 



<DD, 



(D*Dq), 



(*) 7-)(*)l 



BD ' ,D 



-if D *m D * 



)(m B + m D *)A 1 (m 2 ) 



-{e* D * -p D *)(e* D * ■ p D * ) mB+m q D , 

2V(m 2 „) 



(16) 



For the penguin contributions, we will consider not only QCD and electroweak penguin 
operator contributions but also the contributions from the electromagnetic and chromo- 



magnetic dipole operators and Q$ g , as defined by the factor Pf [22 



Pi = 0, 
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y- G D (*»K) 



(17) 



with the penguin loop-integral function G D (*)(m p ) defined as 



G „(.)(m p ) 
G(m p , k) 



duG(m p , /c)$„(*) (ii), 



-4 / <ixx(l — x)ln 
Jo 



m 2 — k 2 x(l — x) 



mt 



(18) 
(19) 



where k 2 = m 2 c + u(m 2 — m 2 c — m\ l2 ) + u 2 m 2 M2 is the penguin momentum transfer with 
u = 1 — u. In the function G n (*)(m p ), we have used a meson-emitting distribution 
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amplitude <£>(*> (tt) = 6u(l—u)[l+a n (*)(l — 2u)], in stead of keeping k 2 as a free parameter 



as usual. The constants F% and F 2 in Eq. (TT7|) are defined by [22] 

" l —um 2 n —2m 2 +m b m c 



Jo du$ Dq (u 
So du ® 

So du ®D q 



m b —m c 
m b 
~W 



2um 
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j'l du$ Dq (u)^[u(m b - m c ) + m c ] 


Jq 1 du$ Dq {u)^[u(m b + m c 




(DD q ), 
{DDI), 

(D*Dg), 

(D*D*), 



{DD* q ), 

{D*D q ), 
(D*D* q ), 



(20) 



(21) 



where e 2L ■ p\ 



[mi 



m 2 M , - ml) / (2m M *) and e* 2T ■ p x = for B -> D*D* decays. 



C. Observables of B — > M1M2 decays 

In the B meson rest frame, the branching ratios of two-body B meson decays can be 
written as 



B(B^D M D M ) = ^h 



57T m 



B 



M(B ->■ D^D^) 2 , (22) 



where tb is the B meson lifetime, and \p c \ is the magnitude of momentum of particle Mi 
and M 2 in the B rest frame and written as 

J[m 2 B - (m DW +m DW ) 2 ][m| - (m D <.) -m D w) 2 ] 

^ = " -2^ — < 23 > 

In 5 — > D*D* decays, one generally should evaluate three amplitudes as M.o,± m the 
helicity basis or as Ml,\\,± i n the transversity basis, which are related by M. L = M. and 
M\\± = M+± J*- . Then we have 

II'- 1 - V2 

\M(B ->■ D*D*)| 2 = |7W | 2 + l-M+T + l-M-l 2 = |-Ml| 2 + l^nl 2 + \M ± \ 2 . (24) 

The longitudinal polarization fraction fi and transverse polarization fraction f± are 
defined by 

f L , ± (B D*D*) = ^ = , kA ,- '^'if, l2 - (25) 



In charged i? meson decays, where mixing effects are absent, the only possible source 
of CPAs is 

^, dir = \M k (B-^J)/M k (B+^f)\ 2 -l 
° P |M fc (S-^7)/M fc (S+^/)| 2 + l' 



7 



and k 



for B — > D*D* decays and k = L for B u — > DD q ,DD*D*D q decays. 



Then for B u — >■ D*D* decays, we have 



A^f{B -> D*D*) 



12 i /J^,dir| A/( |2 



l-^lll 2 



(27) 



For neutral -B 9 meson decays, the situation becomes complicated because of B® — B® 
mixing, and have been studied by many authors. We do not repeat the lengthy discussions 
here, one can see Refs. HIl for details. 



IV. NUMERICAL CALCULATIONS 
A. Input parameters 

• CKM matrix elements: In numerical calculation, we will use the following values 



which given as [39 



v ud \ 


= 0.9743, 


\v us 


= 0.2252, 


\v ub 


= 0.0035, 


\Vcd\ 


= 0.2251, 


\Vcs\ 


= 0.9735, 


\Vcb\ 


= 0.0412, 


\Vtd\ 


= 0.0086, 


\Vts\ 


= 0.0404, 


\Vtb\ 


= 0.9991, 



P = (21.58 



+0.91\o 
0.8lJ > 
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(67.8j42)°. 



(28) 



Quark masses. When calculating the decay amplitudes, the pole and current quark 
masses will be used. For the former, we will use 



m u = 4.2MeV, m c = 1.5GeV, m t = 175GeV, 
m s = 0.122GeV, m h = 4.62GeV. 



m d = 7.6MeV, 



The current quark mass depends on the renormalization scale. In the MS scheme 
and at a scale of 2GeV, we fix 

m u (2GeV) = 2.4MeV, m d (2GeV) = 6MeV, 

m s (2GeV) = 105MeV, m b (fn b ) = 4.26GeV, 



and then employ the formulae in Ref . |33 
m(/i) =m(fi ) "7 r '- 1 





^ r 


a 8 (fi) 


2/3 


_a s (/i )_ 





,(°) 



2/?o 



2/?o 2 



Ait 



(29) 



to obtain the current quark masses at any scale. The definitions of a s , 7m\ jm , Po, 
and Pi can be found in Ref. (33| . 

Decay constants: The decay constants of D* mesons have not been directly measured 
in experiments so far. In the heavy-quark limit (m c — > oo), spin symmetry predicts 
that fo* = fog-, and most theoretical predictions indicate that symmetry-breaking 
corrections enhance the ratio fu* / fjj by 10% — 20% |40j, |4l|. In this paper, we will 
take f D = 0.201 ± 0.017GeV, fo s = 0.249 ± 0.016GeV and f D * = f Dq as our input 
values. 
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Distribution amplitudes: The distribution amplitudes of mesons are less con- 
strained, and we use the shape parameter a D (*) = 0.7 ± 0.2 and a m = 0.3 ± 0.2. 

D s 

Form factors: For the form factors involving B — > D" transitions, we take expres- 
sions which include perturbative QCD corrections induced by hard gluon vertex cor- 
rections of b — y c transitions and power corrections in orders of l/rrib^ (HI E2}. As 
for Isgur-Wise function we use the fit result £(a;) = 1 — 1.22(u; — l)+0.85(u; — l) 2 
fromRef. 0. 



Mass and lifetimes: For B and D meson masses, the lifetimes, we use the following 
as input parameters 44 . 





= 5.279GeV, 


B d 


= 5.280GeV, 


m Bg = 5.366GeV, 




= 1.865GeV, 


M D+ 


= 1.870GeV, 


M D + = 1.969GeV, 




= 2.007GeV, 




- = 2.010GeV, 


M m+ = 2.107GeV 




= (1.638)ps, 




(1.530)ps, 




= (1-425™ 


)ps. 







(30) 

Using the input parameters given above, we then present the numerical results and 
make some theoretical analysis for double charm B u <i and B s decay processes. 



B. data and theoretical prediction 

1. b — > ccd decays 

In the SM, B° d D M+ D w ~ , B~ -)• D M ° D M ~ and B° s -)• D { s " )+ D M ~ decays are 
dominated by the tree b — > ccd transition, and receive additional b — > ccd penguin diagram 
contributions. 

In Table HTl we show the theoretical predictions for the CP-averaged branching ratios 
and the polarization fractions in SM and mSUGRA model. The weighted averages of 
the relevant experimental data [lij are given in the last column in both the Table [TT] and 
Table IIHI The data with a star in the top right corner denote the BaBar measurement 
only, while that with two stars are the Belle measurements only. The central values of the 
theoretical predictions are obtained at the scale /1 = nib, while the two errors are induced 
by the uncertainties of f D = 0.201 ± 0.017GeV and 7 = 67.8° ± 20°. 

From the numerical results and the data as given in Table [Til we have the following 
remarks on the branching ratios and the polarization fractions of b — > ccd double charm 
decays: 

(i) The SUSY contributions to the branching ratios of the considered decays are 
indeed very small, less than 5%, which is consistent with the general expectation 
since these decays are all "tree" dominated decay processes. 

(ii) Thhe theoretical predictions of the Br's in both the SM and the mSUGRA model 
are consistent with the experimental measurements within ±2<r errors. The central 
value of the theoretical prediction for Br(B® — > D + D~) (Br(B~ — > D*°D~)) is, 
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however, much larger (smaller ) than that of the corresponding measurement. This 
point will be clarified by the forthcoming LHC experiments. 

(iii) The SUSY contributions to the polarization fractions of these decays in 
mSUGRA model are very small, less than 2%, and can be neglected safely. Only 
the central values are presented here since they are not sensitive to the variations 
of the form factors and the weak phase 7, which can be seen from the definition of 
the polarization fraction. 



TABLE II: Theoretical predictions for CP-averaged branching ratios (in units of 1CP 4 ), polar- 
ization (in percent) for b — > ccd decays in the SM and mSUGRA model. The last column shows 
currently available data fi^ . 



Observables 



SM 



mSUGRA 



(A) 



(B) 



Data 



B(S° d 



D+D-) 
£)*±.D=F) 

D* + D*-) 



o 9fi +0.57- 

tr O9+1-05- 
°- yz -0.95- 
7 94+1-28- 
' - z ^-1.17- 



-0.10 
-0.12 
-0.01 
-0.01 
-0.06 
-0.06 



0.58H 
-0.53- 
1.04H 
0.96- 
1.28H 
1.17- 



0.10 
0.11 
0.01 
0.01 
0.06 
0.06 



+0. 
-0. 

+1. 



19±i 



55+0.08 
51-0.13 
04+0.01 
96-0.01 
26+0.06 
17-0.06 



2.1 ±0.3 

6.1 ± 1.5 

8.2 ±0.9 



B{B- 
B(B- 
B(B~ 



D°D~) 
D*°D~) 
D°D*~) 
D*°D*~ 



3A3~ 



3.48j££j 

+0.60H 
0.51- 
O no-l-0.5lH 
z - yz -0.15- 

' • ' °-1.16- 



0.11 
0.78 
0.03 
0.07 
0.02 
0.03 
0.05 
0.07 



0.62H 
0.57- 
0.60H 
0.56- 
O9+0.52H 
yz -0.47- 
1.36-1 



50] 

43 H 



76] 



1.26- 



0.10 
0.12 
0.02 
0.02 
0.02 
0.02 
0.06 
0.07 



37: 

44j 
.89^ 



.681} 



59+0.11 

55-0.14 
61+0.03 
55-0.02 
51+0.03 
47-0.03 
36+0.07 
22-0.07 



3.8 ±0.4 
6.3 ±1.4 ±1.0* 

3.9 ±0.5 
8.1 ± 1.2 ± 1.2* 



B(B° S 
B(B° S 
B(B° S 
B(B° S 



D+D-) 
D*+D~ 
D+D*- 
D*+D*- 



3.22 



+0.51- 



-0.52- 



7.12 



0.10 
0.11 
0.02 
0.02 
0.03 
0.02 
0.07 
-1.15-0.06 



O I q + 0.55 

°- 1J -0.51 

2 67 +0 - 48 
-0.43 



-1.26H 



24 



-0.57-1 
-0.53- 



13 



0.09 
0.12 
0.02 
0.02 
0.02 
0.03 
0.06 
-1.15-0.06 



1 q + 0.55 

1J -0.51 
ro+0.47 



-1.26H 



HIS 

1418 
651° 
+1 



07 



55+0.11 

50- 0.13 
55+0.02 

51- 0.02 
47+0.02 
45-0.03 
24+0.06 
15-0.08 



h(B° d 

Mb- 

h(B° s 



D* + D*-) 
■ D*°D*-) 
D*+D*-) 



53.86 
53.88 
53.88 



53.87 
53.89 
53.89 



53.79 
53.81 
53.81 



57.0 ±8.0 ±2.0* 



h(B° d 
f±(B- 

f±(B°s 



D*+D*- 
► D*°D*- 
D* + D*- 



5.51 
5.52 
5.20 



5.50 
5.52 
5.20 



5.51 
5.53 
5.21 



15.0 ±2.5 



In Table IIII1 we present the theoretical predictions for the CPAs in the framework of 
the SM and the mSUGRA model. The currently available data are also listed in the last 
column. The uncertainties come from the scale m&/2 < \i < 2m{, and the weak angle 
7 = 67.8° ± 20°. From the numerical results and the data, we find that 

(i) Just as generally expected based on the SM, the direct CPAs Cf are indeed quite 
small, while the mixing-induced CPAs of all considered decays are close to —0.7: 
i.e. S f w sin(2/3) w -0.7. 
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TABLE III: Theoretical predictions of CPAs (in percent) for the exclusive color-allowed b 
decays. The last column shows the word averages 44]. 



ccd 



Observables 



SM 



mSUGRA 



(A) 



(B) 



Data 



D+D-) 
D*+D~) 
D+D*-) 
-> D* + D*- 



7 c Q+l.4+1.4 
"' J,o -1.5-0.6 



-68.4" 



-0.2+0.3 
-0.3-0.2 



-75.li^ 



-68.4" 



RS /|+0. 1+0.2 
-OO-4-0.4-0.2 

7n 9 +0.4+0.4 
~ ' u -^-0. 6-0.1 



-1.3 
-0.6 
,+0.2+0.3 
-0.3-0.2 
RS 4+0.1+0.2 
-t)8.4_ 04 _ 02 

7n 1+0.5+0.3 
" ' u - 1 -0.6-0.2 



-76.3ijj 



-68.5 
-68.5 



-70 4+0- 4 +0- 4 

' u - 4 -0.7-0.2 



1.6 
0.7 
+0.2+0.3 
0.3-0.2 
+0.2+0.2 
0.4-0.2 



-87 ± 26 
-61 ± 19 
-78 ±21 
-81 ± 14 



C + (B° d ,B° d 



D+D-) 
D*+D~) 
D+D*-) 
-> D* + D*- 



-4.4 
7. 

-8.4 
-1.2 



+0.3+1.0 
0.4-0.5 
+0.3+0.7 
-0.6-0.6 
1.1+0.7 
1.1-0.8 
+0.2+0.2 
0.4-0.1 



+0.3+1.0 
-0.4-0.5 
7 7+0.3+0.7 
'•'-0.6-0.6 
j+1.1+0.7 
-1.1-0.8 
+0.2+0.2 
-0.4-0.1 



-4.4 



-8.3] 
-1.2 



+0.3+1.0 
-0.4-0.6 
o q+0.3+0.8 
°-°-0.6-0.7 
-,+1.0+0.8 
-1.0-0.9 
,+0.2+0.2 
-0.4-0.1 



-4.5 



-8.9: 

-1.2^ 



-48 ± 42 
-9 ±22 
7 ±14 
-7 ±9 



■Acp{B u 
^cp(B u ■ 
■A-q?{B u ■ 
A+' dir (B- 



D°D~) 
D*°D-) 
D°D*~) 
D*°D*- 



a 4+0.4+1.0 
^-0.3-0. 2 



-0.6 
1.2 



+0.4+0.1 
0.2-0.1 
+0.4+0.1 
0.2-0.2 



1.2 



+0.4+0.1 
-0.2-0.2 



4 4+O.4+O 
^•+-0.3-1 

-0.6 
1.2 
1.2 



5 


0.4+0.1 
0.2-0.1 
+0.4+0.1 
0.2-0.2 
+0.4+0.1 
0.2-0.2 



4 rr+0. 4+0.6 
^•°-0. 3-1.0 



-0.6 
1.2 
1.2 



-0.4+0.1 
-0.2-0.1 
+0.4+0.1 
0.2-0.2 
+0.4+0.1 
-0.2-0.2 



-3 ±7 
13 ± 18 ± 4* 

3± 10 
-15 ±11 ±2* 



A^(B° S - 

^M4° - 

^ t p(B° s - 
^P dir (4° 



D+D") 

D*+D~) 

D+D*') 



4.4 
-0.6 



+0.4+0.5 
0.3-1.0 
+0.4+0.1 
0.2-0.2 



4.4 
-0.6 



+0.4+0.5 
0.3-0.6 
+0.4+0.1 
0.2-0.1 



4.5 
-0.6 



+0.4+0.6 
0.3-1.0 
+0.4+0.1 
0.2-0.1 



I 2+0.4+0.1 



1.2 



-0.2- 
-0.4- 



0.2 
0.1 



-0.2-0.2 



i 9+0.4+0.1 
-■-■^-0.2-0.2 
-, 9+0.4+0.1 
-■-• z -0.2-0.2 



1 9+0.4+0.1 
-■■■ z -0. 2-0.2 



1.2 



+0.4+0.1 
-0.2-0.2 



(ii) The SUSY contributions to all considered decays are less than 7%. The new 
physics contributions is not sensitive to the variation of the scale \i and the weak 
angle 7. 

(iii) The theoretical predictions in the SM and mSUGRA model are all consistent 
with the experimental measurements within ±la error. Of course, the errors of 
currently available data are very large now. 



2. b — > ccs decays 

The twelves decay modes — > D ( '* )+ D { * ] , B~ — > D { * >0 and B Q S — > D ( * )+ are 
the tree-dominated processes, and also receive the additional b — > ccs penguin contribu- 
tions. 

In Table HV] we present the theoretical predictions for the CP-averaged branching ratios 
and the polarization fractions in the framework of the SM and the mSUGRA model. The 
last column in table HVl correspond to the world averages [44j ]. The theoretical predictions 
for CP asymmetries of considered decays are given in Table |Vl although they have not been 
measured yet. The central values of the theoretical predictions are obtained at the scale 
/1 = m,;,, while the two errors are induced by the uncertainties of fo = 0.201 ± 0.017GeV 
and 7 = 67.8° ± 20°. 
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TABLE IV: Theoretical predictions for CP-averaged £>V (in units of 10 -3 ) and polarization 
fractions (in units of 10 -2 ) of exclusive color-allowed b — > ccs decays in the SM and the mSUGRA 
model. The last column corresponds to the world averages 44]. 



Observables 



SM 



mSUGRA 



(A) 



(B) 



Data 



B(S° d 
B(B° d 



D+D-) 
D*+D-) 

D+D*-) 
D*+D*-) 



8.771 1 



8.78 
7.30 
21.2 



4-1 
1 



+2. 



16+0.02 

09- 0.02 
16+0.01 

10- 0.01 
97+0.01 
91-0.01 

I ±0.0 



+1 
1 
+1 

7.3llg 

>+2 



3.83 
3.77 



17+0.02 
10-0.02 
17+0.01 
09-0.01 
97+0.01 
91-0.01 
8 



8.39 
8.78 
7.22 



1.11+0.02 
1.36-0.02 
+1.17+0.01 
1.09-0.01 
0.96+0.01 
0.90-0.01 



21.2 + S« ±0.0 20.9 



+2.8 
-2.6 



±0.0 



7.4 ±0.7 
8.2 ± 1.1 

7.5 ± 1.6 
17.8 ± 1.4 



B(B- 
B{B~ 



D°D 7 ) 
D*°D-) 
D°D*~) 
D*°D*~) 



9.38 



+i 
l 

9.401J 
7.82+J 
22.6l| 



24+0.01 
17-0.02 
24+0.01 
17-0.01 
04+0.01 
97-0.01 

8 + 0-0 



9.44 



+i 
l 

9.391} 
7.83l^ 
22.7l| 



25+0.01 
18-0.02 
25+0.01 
17-0.01 
04+0.01 
97-0.01 

8 + 0-0 



8.97 



+i 
l 

9.401J 
7.731J 
22.41 3 ; 



.19 
.12 
.25 
.17 
.0.3 
.96 

°±0.0 



^0.02 
-0.02 
1-0.01 
-0.01 
1-0.01 
-0.01 



10.2 ± 1.7 
8.4 ± 1.7 
7.8 ± 1.6 
17.4 ±2.3 



B(B°s 
B(B° S 

B(B° S 



DfD-) 
D* s + D-) 
DfD* s -) 
D*+D*-) 



1.68 



8.741* 
7.16+° 
20.8 



+2 



15+0.02 

08- 0.02 
16+0.01 

09- 0.01 
95+0.01 

0.01 

8 



3.73 



-0.89 



±0.0 



3.731* 

+o 
-o 

+2. 



7.17 
20.8 



16+0.02 
08-0.02 
16+0.01 
08-0.01 
98+0.01 
88-0.01 
8 



8.30 
8.75 
7.08 
±0.0 20.6 



10+0.02 
03-0.02 
16+0.01 
09-0.01 
94+0.01 
88-0.01 



+2.7 
6 



±0.0 



11 ±4 

< 121 

< 257 



/urn 

h(B- 

h(B° s 



D* + D*~) 



D*°D 



s 

D* s + D*-) 



51.68 
51.70 
51.70 



51.70 
51.72 
51.71 



51.58 
51.61 
51.60 



52 ±5 



f±(B° d 
f±(B~ 

f±(B°s 



D* + D*~) 
■ D*°D*-) 
D*+D*~) 



5.50 
5.51 
5.19 



5.50 
5.51 
5.18 



5.51 
5.52 
5.20 



From the numerical results and currently available data, one can see that 

(i) For the Br's and CPAs, the SUSY contributions again are very small for all 
considered decays, less than 3% numerically. The theoretical predictions in both 
the SM and the mSUGRA model are all consistent with currently available data 
within one or two standard deviations. 

(ii) The direct CP violations C{B° S -»■ D* + D-) and C(B° S D+D*~) are at the 
±10% level and to be tested by the LHC experiments. And the CP asymmetries 
for the remaining ten decays are very small, about 10~ 3 or 10~ 4 numerically, since 
the penguin effects are doubly Cabibbo-suppressed for the color-allowed b — >■ ccs 
decays. 



V. SUMMARY 

In this paper, we have investigated the new contributions to the branching rations, 
polarization fractions and CP asymmetries of the twenty three double charm decays 
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TABLE V: Theoretical predictions for CPAs (in percent) of exclusive color-allowed b — > ccs 
decays in the SM and the mSUGRA model. 



Observables 



SM 



mSUGRA 



(A) 



(B) 



Data 



i dir (R° 



D+D-) 

D*+D~) 

D+D*-) 



_fl or+O.02+0.05 
u - zu -0.03-0.02 

n n o+o.02+o.oi 

u - uo -0.02-0.01 
+0.02+0.02 
0.02-0.01 



-0.07 



_n n7+°- 02 + - 02 

u - u ' -0.02-0.01 



_fl OR+0.02+0.05 
u - zu -0.01-0.02 

n n o+o.02+o.oi 

u - UJ -0.02-0.01 

_n n7+°- 02 + - 02 
u - u ' -0.01-0.01 

_n n7+°- 02 + - 02 
u - u ' -0.01-0.01 



_n O7+0.02+0.06 
u - z ' -0.02-0.02 

n n o+o.02+o.oi 

u - uo -0.02-0.01 

_n n7+°- 02 + - 02 

u - u ' -0.02-0.01 

_n n7+°- 02 + - 02 

u - u ' -0.02-0.01 



Acp{B u 

A(yp(B u 

Ao?{B~ 

x+,dir/ R - 



D*°D-) 
D°D*~) 
+ D*°D*-) 



_fl 9R+0.02+0.05 
u - zu -0.03-0.02 
q Q3+0. 02+0.01 



-0.26 



+0.02+0.05 
-0.01-0.02 



-0.27 



+0.02+0.06 
-0.02-0.02 



-0.02-0.01 
+0.02+0.02 
-0.02-0.01 
n n7+0-02+0.02 
-U.Urf _q. 02-0. 01 



n n o+o.02+o.oi 

u - uo -0.02-0.01 



n n o+o.02+o.oi 

u - uo -0.02-0.01 



-0.07 



-0.07 



+0.02+0.02 
-0.01-0.01 
n n7 +0.02+0.02 
" U - U( -0.01-0.01 



-0.07 



-0.02+0.02 
-0.02-0.01 
n n7 +0.02+0.02 
~ u - u ' -0.02-0.01 



S(B° S ,B° S ^ 
S + (B° S ,B° S - 



D+D-) 
D*+D-) 
D+DT) 
> D*s + D*-) 



p, Cq + 0.11 + 0.11 

U -' JO -0.12-0.12 
fl nq+0.02+0.02 
u - yo -0.01-0.02 
n Q ,+0.03+0.02 
u - y4 -0.01-0.01 

n i q+o.04+0.01 

u - io -0.04-0.03 



+0.06+0.04 
0.11-0.10 
+0.01+0.02 
0.06-0.02 
n QA +0. 10+0.02 
_u - y4 -0.01-0.02 



0.51 
0.93 



n R9 +0.11+0.06 
u -°^-0.11-0.12 
n Q/i+0.02+0.02 
u - y ^-0.01-0.02 
n nq+0.02+0.02 
u - yj -0.03-0.02 



0.12 



+0.03+0.01 
-0.03-0.02 



0.14 



+0.05+0.02 
-0.03-0.02 



C(B° a ,B° a - 
C(B° a ,B° a - 
C(B° S ,B° S - 

c+(Bim 



DfD-) 

DfD*-) 
+ D*+D* s 



0.26 
9.91 



+0.03+0.05 
-0.02-0.02 
+0.91+0.05 



-9.93 
0.07 



-1.14-0.04 
+1.16+0.01 
-0.95-0.04 
+0.01+0.01 
-0.02-0.02 



0.26 
9.82 



-0.01+0.02 
-0.02-0.05 
+0.21+0.04 
-1.15-0.05 



0.27 



+0.02+0.02 
-0.02-0.06 



10.52 



+0.89+0.05 
-1.12-0.05 



_q 84 +l. 18+0.05 
y -°^-0.25-0.03 



■10.54 



+1.14+0.05 
-0.93-0.05 



0.07 



+0.02+0.01 
-0.01-0.02 



0.07 



+0.01+0.01 
-0.02-0.02 



B/B s — >■ DyIdyI in the SM and the mSUGRA model by employing the effective hamil- 
tonian for AB = 1 transition and the naive factorization approach. 

From the numerical results and the phenomenological analysis, the following conclu- 
sions can be reached: 

(i) For the exclusive double charm decays B/B s — > D^D^ studied in this pa- 
per, the SUSY contributions in the mSUGRA model are very small, less than 7% 
numerically. It may be difficult to observe so small SUSY contributions even at 
LHC. 

(ii) All the theoretical predictions in the SM and mSUGRA model are consistent 
with the experimental measurements within ±2cx errors. 

(iii) The theoretical predictions in both the SM and mSUGRA model still have large 
theoretical uncertainties. The dominant errors are induced by the uncertainties of 
the form factors fp or f Ds . 



Acknowledgments 

We are grateful to Wen-juan Zou for valuable help. This work is partially supported 
by the National Natural Science Foundation of China under Grant No. 10947020, and by 



13 



Foundation of Henan Educational Committee for Youth Backbone Scholars in Colleges 
and Universities, and by the Natural Science Foundation of the Eduction Department of 
Henan Province under Grant No. 2010A140012. 



S. Fratina et al, Phys. Rev. Lett. 98, 221802 (2007). 

J. Anderson et al. [BABAR Collaboration], Phys. Rev. D 79, 032002 (2009). 
T. Aushev et al. [Belle Collaboration], Phys. Rev. Lett. 93, 201802 (2004). 
B. Aubert et al. [BABAR Collaboration], Phys. Rev. D 76, 111102 (2007). 
H. Miyake et al. [Belle Collaboration], Phys. Lett. B 618, 34 (2005). 
K. Vervink, larXiv:0810.3l67l [hep-ex]. 

B. Aubert et al. [BABAR Collaboration], Phys. Rev. D 73, 112004 (2006). 
K. Abe et al. [Belle Collaboration], |ar2Civi0708 i 1668J [hep-ex]. 

G. Majumder et al. [Belle Collaboration], Phys. Rev. Lett. 95, 041803 (2005). 
Run-Hui Li, Xiao-Xia Wang, A.I. Sanda and Cai-Dian Lii, Phys. Rev. D 81, 034006 (2010). 
Ying Li, Cai-Dian Lii and Zhen-Jun Xiao, J. Phys. G 31, 273 (2005). 
R. Zwicky, Phys. Rev. D 77, 036004 (2008). 
R. Fleischer, Eur. Phys. J. C 51, 849 (2007). 

M. Gronau, J. L. Rosner and D. Pirjol, Phys. Rev. D 78, 033011 (2008). 

C. S.Kim, Ru-min Wang and Ya-Dong Yang, Phys. Rev. D 79, 055004 (2009). 

H. P. Nills, Phys. Rep. 110, 1 (1984), and references therein. 

H. Baer, M. Brhlik, D. Castano and X. Tata, Phys. Rev. D 58, 015007 (1998); H. Baer 
and M. Brhlik, Phys. Rev. D 55, 3201 (1997). 

T. Goto, T. Nihei and Y. Okada, Phys. Rev. D 53, 5233 (1996); ibid 54, 5904 (Erratum) 
(1996). 

T. Goto, Y. Okada, Y. Shimizu and M. Tanaka, Phys. Rev. D 55, 4273 (1997); ibid 66, 
019901 (Erratum) (2002). 

C.S. Huang and X.H. Wu, Nucl. Phys. B 657, 304 (2003). 

Z.J. Xiao and W.J. Zou, Phys. Rev. D 70, 094008 (2004); W.J. Zou and Z.J. Xiao, Phys. 
Rev. D 72, 094026 (2005); W.J. Zou and Z.J. Xiao, Phys. Rev. D 73, 057504 (2006). 
M. Wirbel, B. Stech and M. Bauer, Z. Phys. C 29, 637 (1985); M. Bauer, B. Stech and 
M. Wirbel, Z. Phys. C 34, 103 (1987). 

A. Ali, G. Kramer, and CD. Lii, Phys. Rev. D 58, 094009 (1998); HY. Cheng, Phys. 
Lett. B 335, 428 (1994); HY. Cheng, Z. Phys. C 69, 647 (1996). 

M. Beneke, G. Buchalla, M. Neubert, and C.T. Sachrajda, Phys. Rev. Lett. 83, 1914 
(1999); Nucl. Phys. B 591, 313 (2000). 

M. Beneke and M. Neubert, Nucl. Phys. B 651, 225 (2003). 

H.-n. Li and H.L. Yu, Phys. Rev. Lett. 74, 4388 (1995); Phys. Lett. B 353, 301 (1995); 
H.-n. Li, Phys. Lett. B 348, 597 (1995); H.-n. Li and H.L. Yu, Phys. Rev. D 53, 2480 
(1996). 

[27] R. Barbieri, S. Ferrara, and C.A. Savoy, Phys. Lett. B 119, 343 (1982); A.H. Chamseddine, 
R. Arnowitt, and P. Nath, Phys. Rev. Lett. 49, 970 (1982); L. Hall, J. Lykken, and 
S. Weinberg, Phys. Rev. D 27, 2359 (1983); Nobuyoshi Ohta, Prog. Theor. Phys. 70, 542 



(1983); A. Djouadi, et al, hep-ph/9901246 



14 



Gi-Chol Cho and K. Hagiwara, Nucl. Phys. B 574, 623 (2000). 

J. Ellis, S. Heinemeyer, K.A Olive, and G. Weiglein, JHEP 0502, 013 (2005); S. Heinemeyer, 
W.Hollik and G. Weiglein, |hep-ph/041221~4l 



For more details of SPA project, see SPA home page: http://spa.desy.de/spa/ 
B.C. Allanach, et al, Eur. Phys. J. C 25, 113 (2002). 

A. Ali, E. Lunghi, C. Greub, and G. Hiller, Phys. Rev. D 66, 034002 (2002). 

G. Buchalla, A. J. Buras and M. E. Lautenbacher, Rev. Mod. Phys. 68, 1125 (1996). 

M. Neubert and V. Rieckert, Nucl. Phys. B 382, 97 (1992). 

M. Gronau, Phys. Lett. B 233, 479 (1989). 

J. Soto, Nucl. Phys. B 316, 141 (1989). 

W. F. Palmer and Y. L. Wu, Phys. Lett. B 350, 245 (1995). 

A. Ali, G. Kramer and C. D. Lii, Phys. Rev. D 59, 014005 (1998). 

J. Charles et a/.(CKMfitter Group), Eur. Phys. J. C 41, 1 (2005), updated results and plots 



available at: http://ckmfitter.in2p3.fr/ 



M. Neubert, Int. J. Mod. Phys. A 11, 4173 (1996). 
M. Neubert, Phys. Rept. 245, 259 (1994). 
M. Neubert, Phys. Rev. D 46, 2212 (1992). 

H. Y. Cheng, C. K. Chua and C. W. Hwang, Phys. Rev. D 69, 074025 (2004). 
C. Amsler et al, (Particle Data Group), Phys. Lett. B 667, 1 (2008). 



15 



